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Abstract Block elimination algorithms for solving sparse discrete optimiza- 
tion problems are considered. The numerical example is provided. The bench- 
marking is done in order to define real computational capabilities of block 
elimination algorithms combined with SYMPHONY solver. Analysis of the 
results show that for sufficiently large number of blocks and small enough 
size of separators between the blocks for staircase integer linear program- 
ming problem the local elimination algorithms in combination with a solver 
for solving subproblems in blocks allow to solve such problems much faster 
than used solver itself for solving the whole problem. Also the capabilities of 
postoptimal analysis (warm starting) are considered for solving packages of 
integer linear programming problems for corresponding blocks. 



1 Introduction 

The use of discrete optimization (DO) models and algorithms makes it possi- 
ble to solve many practical problems, since the discrete optimization models 
correctly represent the nonlinear dependence, indivisibility of an objects, con- 
sider the limitations of logical type and all sorts of technology requirements, 
including those that have qualitative character. But unfortunately, most of 
the interesting problems are in the complexity class TVP-hard and may require 



A. Sviridenko 

Faculty of Mathematics and Computer Science 

Tavrian National University, Vernadsky Av. 4, Simferopol 95007 Ukraine 
e-mail: |oleks . svi ridenko@gmail . com| 

O. Shcherbina 

Faculty of Mathematics, University of Vienna 
Nordbergstrasse 15, A-1090 Vienna, Austria 
e-mail: |oleg . shcherbinaBunivie . ac . at | 



1 



2 



Alexander Sviridenko and Oleg Shcherbina 



searching a tree of exponential size in the worst case. Many real DOPs from 
OR applications contain a huge number of variables and/or constraints that 
make the models intractable for currently available solvers. Usually, DOPs 
from applications have a special structure, and the matrices of constraints for 
large-scale problems have a lot of zero elements (sparse matrices), and the 
nonzero elements of the matrix often fall into a limited number of blocks. The 
block form of many DO problems is usually caused by the weak connectedness 
of subsystems of real- world systems. 

Among the block structures let us pay particular attention to the block- 
tree structure, a special case of which is a staircase or quasiblock structure 
(Fig. [J). The problems of optimal reservation of hotel rooms [2] have the 
quasiblock structure, similar to the previous temporal knapsack problems |10| . 
recently received the application in solving problems of the prior reservation 
of computing resources in the Grid Computing. 



Fig. 1 Quasi-block structure. 

One of the promising ways to exploit sparsity in the constraint matrix 
of DO problems are local elimination algorithms (LEA)p], including local 
decomposition algorithms nonserial dynamic programming (NSDP) [6] 
algorithms, 0. To extract special block structures there are such promising 
graph-based decomposition approaches as methods of tree decomposition '5^ . 
The purpose of this paper to define real computational capabilities of block 
elimination algorithms combined with modern solvers. 



2 Local elimination algorithms for solving discrete 
problems 

2.1 General scheme of local elimination algorithms 

In the papers [9], [5] considered the general class of local elimination algo- 
rithms for computing information, that have decomposition approach and 
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that allow to calculate some global information about a solution of the entire 
problem using local computations. 

A local elimination algorithm (LEA) eliminates local elements of the prob- 
lems structure defined by the structural graph by computing and storing local 
information about these elements in the form of new dependencies added to 
the problem. 

The local elimination procedure consists of two parts: 

1. The forward part eliminates elements, computes and stores local solutions, 
and finally computes the value of the objective function; 

2. The backward part finds the global solution of the whole problem using 
the tables of local solutions; the global solution gives the optimal value 
of the objective function found while performing the forward part of the 
procedure. 

The algorithmic scheme of the LEA is a directed acyclic graph (DAG) in 
which the vertices correspond to the local subproblems and the edges reflect 
the informational dependence of the subproblems on each other. 

It is important that aforementioned methods use just the local informa- 
tion (i.e., information about elements of given elements neighborhood) in a 
process of solving discrete problems. Thus local elimination algorithms allow 
to calculate some global information about a solution of the entire problem 
using local computations. 

The structure of discrete optimization problems is determined either by the 
original elements (e.g., variables) with a system of neighborhoods specified for 
them with help of structural graph and with the order of searching through 
those elements using a LEA or by various derived structures (e.g., block or 
tree- block structures). 



2.2 Local elimination algorithms 

Consider LEA in details for solving sparse problems of integer linear pro- 
gramming in the case when structural graph is interaction graph of variables, 
which is also called constraint graph. 

Consider the integer linear programming (ILP) problem with binary vari- 
ables 



n 




(1) 



subject to constraints 



n 




(2) 



4 



Alexander Sviridenko and Oleg Shcherbina 



X, =0,l,J = l,2,...,n. (3) 

Definition 1. Variables x and y interact in ILP problem with constraints if 
they both appear in the same constraint. 

Definition 2. Interaction graph of the ILP problem is an undirected graph 
G = iX,E), such that 

1. Vertices X of G correspond to variables of the ILP problem; 

2. Two vertices of G are adjacent iff corresponding variables interact. 

Further, we shall use the notion of vertices that correspond one-to-one to 
variables. 

Example 1. Consider an ILP problem with binary variables: 

2xi + 'ix2 -I- a;3 -f 5x4 + 4x5 + Gxg -I- xy — ?• max (4) 

subject to constraints 

3xi + 4x2 + < 6, (5) 

2x2 + 3x3 + 3x4 < 5, (6) 

2x2 +3x5 <4, (7) 

2x3 + 3x6 + 2x7 < 5, (8) 

X, =0,1, J = l,...,7. (9) 

Definition 3. Two vertices x and y oi G — {X,E) are called neighbors if 
(x,y) G E. 

Definition 4. Set of variables interacting with a variable x £ X \s denoted 
by Nb{x) and called the neighborhood of the variable x. For corresponding 
vertices a neighborhood of a vertex x is a set of vertices of interaction graph 
G = {X, E) that are linked by edges with x. Denote the latter neighborhood 
as Nba{x). 

The solution of a sparse discrete optimization problem (1) - (3) whose 
structure is described by an undirected interaction graph G — {X, E) 
with help of LEA was described in pT] in details. Given an ordering 
Xi, X2, . . . , x„, the LEA proceeds in the following way: it subsequently elim- 
inates xi,X2, . . . ,x„ in the current graph and computes an associated local 
information about vertices from hi{Nb{xi)). This process creates a sequence 
of elimination graphs: G° = G,G^, . . . ,G^ , . . . ,G'', where G" = 0. 

The process of interaction graph transformation corresponding to the LEA 
scheme is known as elimination game which was first introduced by Parter 
as a graph analogy of Gaussian elimination. 
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2.3 Block local elimination algorithms 

The local elimination procedure can be applied to elimination of not only 
separate variables but also to sets of variables and can use the so called 
elimination of variables in blocks, which allows to eliminate several variables 
in block. 

Applying the method of merging variables into meta-variables allows to 
obtain condensed or meta-DOPs which have a simpler structure. If the re- 
sulting meta-DOP has a nice structure (e.g., a tree structure) then it can be 
solved efficiently. 

An ordered partition of a set X is a decomposition of X into ordered 
sequence of pairwise disjoint nonempty subsets whose union is all of X. 

In general, graph partitioning is iVP-hard. Since graph partitioning is diffi- 
cult in general, there is a need for approximation algorithms. A popular algo- 
rithm in this respect is MeTi^, which has a good implementation available 
in the public domain. 

An important special case of partitions are so-called blocks. Two variables 
are indistinguishable if they have the same closed neighborhood. A block is a 
maximal set of indistinguishable vertices. The blocks of G partition X since 
indistinguishability is an equivalence relation defined on the original vertices. 
The corresponding graph is called condensed graphs which is a merged form 
of original graph. Not formally, a condensed graph is formed by merging all 
vertices with the same neighborhoods into a single meta-node (supervariable) . 

An equivalence relation on a set induces a partition on it, and also any 
partition induces an equivalence relation. Given a graph G — {X,E), let 
X be a partition on the vertex set X. X — xi,X2, ■ ■ ■ ,Xp, p < n, where 
Xi = Xki{Ki is a set of indices corresponding to xi, I — 1, . . . ,p). For this 
ordered partition X, the DOP (1) - (3) can be solved by the LEA using 
quotient interaction graph G. 

That is, Uf^j^Xi = X and Xj n Xk = for i 7^ k. We define the quotient 
graph of G with respect to the partition X to be the graph 

G = G/X = (X,£), 
where (xj, Xk) £ f if and only if Nbc{:>i-i) fl Xk 7^ 0. 

The quotient graph G(X,f) is an equivalent representation of the inter- 
action graph G{X, E), where X is a set of blocks (or indistinguishable sets of 
vertices) , and f C X x X be the edges defined on X. A local block elimination 
scheme is one in which the vertices of each block are eliminated contiguously. 
As an application of a clustering technique we consider below a block local 
elimination procedure where the elimination of the block (i.e., a subset of 
variables) can be seen as the merging of its variables into a meta- variable. 

A. Forward part 



^ http://www-users.cs.umn.edu/ karypis/metis 
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Consider first the block Xi. Then 

m&x{CNXN\Ais,Xs, <bi, ie M, xj = 0, 1, j € N} = 

max {CN-K^XN-K, + hi{Nb{XK,)\Ais,Xs, <bi, ieM- Ui, 

X, =0,1, j eN-Ki} 
where Ui = {i : Si n Ki ^ 0} and 

hi{Nb{XKj) = max{CK,XKMiS,Xs, <bi,ie Ui, xj = 0, 1, xj G iV6[xi]}. 

The first step of the local block elimination procedure consists of solving, 
using complete enumeration of Xk^ , the following optimization problem 

hi{Nb{XK,)) = ma^{CK^XKMiSiXs, <bi,i€ Ui, xj = 0, 1, xj G Af6[xi]}, 

(10) 

and storing the optimal local solutions Xk^ as a function of the neighborhood 
ofXK„i.e.,X*j,^{Nb{XKj). 

The maximization of f(X) over all feasible assignments Nb{XKi), is called 
the elimination of the block (or meta- variable) Xk^ ■ The optimization prob- 
lem left after the elimination of Xki is: 

max {CN-K,XN-K,+hiiNbiXKj)\Ais,Xs, <bi,iG M-Ui,Xj =0,1, j eN 

X—Xki 

Note that it has the same form as the original problem, and the tabular 
function hi{Nb{XK^)) may be considered as a new component of the modified 
objective function. Subsequently, the same procedure may be applied to the 
elimination of the blocks - meta- variables X2 = Xk2 j • • • : Xp = Xk^ , in turn. 
At each step j the new component /ixj and optimal local solutions X'^ , are 
stored as functions of Nb{XKj \ Xki, ■ ■ ■ ,XKj_i), i-e., the set of variables 
interacting with at least one variable of Xk^ in the current problem, obtained 
from the original problem by the elimination of Xjd, ■ ■ ■ ,XKj_i - Since the 
set Nb{XKp I Xki, ■ ■ ■ ,Xkp_i) is empty, the elimination of Xk,, yields the 
optimal value of objective f{X). 

B. Backward part. 

This part of the procedure consists of the consecutive choice of X'^^, 
-^K ■ ■ ■ i-^Ki^ i-®-' the optimal local solutions from the stored tables 

X*^[~iNb{XK,)),X*^^iNbiXK,\XK,)),...,X*^JXK,_„...,XK,. 

Underlying DAG of the local block elimination procedure contains nodes 
corresponding to computing of functions /ixi (-^6^,(4-1) (xi)) and is a general- 
ized elimination tree. 

Example 2. Consider a DO problem (4) - (9) from example 1 and an or- 
dered partition of the variables of the set X into blocks: Xi = {x5},X2 = 
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{xi, a;2, a;4}, X3 = {a;6,a;7},X4 = {x'a}. For the ordered partition {xi, 2:2, ^3, a;4}, 
this DO problem may be solved by the LEA. Initial interaction graph with 
partition presented by dashed lines is shown in Fig. [5] (a), quotient inter- 
action graph is in Fig. [5] (b), and the DAG of the block local elimination 
computational procedure is shown in Fig. [3] 




a b 

Fig. 2 Interaction graph of the DOP with partition (dashed) (a) and quotient inter- 
action graph (b) (example 2). 




Fig. 3 The DAG (generalized elimination tree) of the local block elimination com- 
putational procedure for the DO problem (example 2). 

A. Forward part 

Consider first the block xi = {X5}. Then Nb{xi) = {2:2}. Solve the fol- 
lowing problem containing 0:5 in the objective and the constraints: 

/ixi (iVfe(xi)) = max{4a;5 | 2x2 + 3x5 < 4, xj £ {0, 1}} 
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and store the optimal local solutions Xi as a function of a neighbor- 
hood, i.e., Xl* (iV6(xi)). Eliminate the block xi. and consider the block 
X2 = {xi,X2,X4}. Nb{x2) = {xs}. Now the problem to be solved is 



subject to 



/ix2 {X3) = max {/ixi {X2) + 2xi + 3x2 + ^Xi} 

3xi + 4x2 + X3 < 6, 

2x2 + 3x3 + 3x4 < 5, 

Xj =0,1, j = 1,2,3,4. 



Build the corresponding table 2. 

<.xi ix2) 



Table 1. 
Calculation of 



Table 2. 

Calculation of h^^{x^) 



X2 


/ixi(a;2) 


2^5 





4 


1 


1 











h^^ix'i) 


xl 


X2 







11 


1 





1 


1 


6 


1 









Eliminate the block X2 and consider the block X3 — {xq^x^}. The neighbor 
of X3 is X3: iV6(x3) = {xa}. Solve the DOP containing X3: 



(2^3) = max{/ix2 + + Sxg + xy | 2x3 + Sxg + 2x7 < 5, Xj £ {0,1}} 

xe,X7 

and build the table 3. 

Table 3. 

Calculation of /ix3(x3) 



X3 




2^6 







18 


1 


1 


1 


12 


1 






Eliminate the block X3 and consider the block X4 = {X3}. Nb{x4) = 0. Solve 
the DOP: 

/ix4 = max{/ix3(a;3),Xj e {0,1}} = 18, 

X3 

where X3 = 0. 

B. Backward part. 

Consecutively find X3*,X2*,Xi*, i.e., the optimal local solutions from the 
stored tables 3, 2, 1. X3 = =4> xg = 1, x^ = 1 (table 3); Xg = ^ xJ = 
1, x^ = 0, xJ = 1 (table 2); x^ = ^ x^ = 1 (table 1). We found the 
optimal global solution to be (1, 0, 0, 1, 1, 1, 1), the maximum objective 
value is 18. 
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3 Research the computational capabihty of local 
algorithms 

3.1 Comparative computational experiment 

Among extremely important research questions about the effectiveness of lo- 
cal elimination algorithms (LEA), the next one causes special interest: "Is 
the use of LEA in combination with a discrete optimization (DO) algorithm 
(for solving problems in the blocks) consistently more efficient than the stan- 
dalone use of the DO algorithm?" [3]. 

Along with the theoretical analysis of performance evaluations, it is of in- 
terest to provide a comparison of LEA combined with other DO algorithms 
by using computational experiments. Providing an exhaustive computational 
study for all possible combinations of LEA with all existing DO algorithms 
(or at least the most efficient ones) is extremely laborious. This work presents 
computational comparisons of the LEA combined with two DO algorithms: 
a) one of the least efficient DO algorithms, which is a simple implicit enumer- 
ation algorithm without use of linear relaxation; b) one of the most effective 
DO algorithm, as the one adopted by the simplex method in the unimodular 
case. 

The purpose of these comparisons is to evaluate the behavior of the LEA 
in combination with very effective algorithms ("lower bound") and weakly 
effective ones ("upper bound"). 

By combining the LEA with implicit enumeration solvers (one of the least 
efficient algorithms) [2] we concluded that for sufficiently large number of 
blocks and small enough size of separators between the blocks for staircase 
ILP problems the performance is better than the stand alone solver. It has 
to be noted, that there are cases that the combination of LEA with the 
simplex method for a small number of blocks performs worse than the simplex 
metho41. Thus, the practical use of LEA in combination with a DO solver 
requires a preliminary machine experiment. Based on this experiment have 
to be determined the acceptable parameters for DO problems. 

The paper [1] presents computational experiment on the effectiveness (in 
terms of precision and time to find a solution) of LEA in combination with 
approximate algorithms from the package of applied programs (PAP) "DIS- 
PRO" developed at the Institute of Cybernetics of Ukrainian Academy of 
Sciences. The studied methods include lexicographic search, recession vector 
and random search. The object of the experiments was a special class of DO 
quasi-block problems, the optimal reservation class, which under certain con- 
ditions are unimodular [2). The results showed that if the efficiency of the 
PAP "DISPRO" algorithms is almost independent of the matrix condition 
structure of the DO problem, but is essentially determined by the dimension 



^ The benchmarking was provided with help of V.V.Matveev. 
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of the problem, then the use of LEA is appropriate for relatively small values 
of the separators and sufficient sparsity of the condition matrix. 

Additionally, by using the LEA in combination with some DO algorithm 
we achieved better accuracy than by using just the DO algorithm. In some 
cases the time to solution is several times lower than the independent calcu- 
lation with an appropriate DO algorithm. By increasing the dimension of the 
problem and by providing the resource time restriction the efficiency of the 
LEA increases. 



3.2 SYMPHONY as a framework for solving mixed 
integer programming problems 

The computational capabilities of the LEA in combination with a modern 
solver were tested by using SYMPHON^vH as the implementation framework. 
SYMPHONY is part of the COIN-OE0 project and it can solve mixed-integer 
linear programs (MILP) sequentially or in parallel. We chose this framework 
since it is open-source and supports warm restarts, which implement postop- 
timal analysis (PA) of ILP problems. 

It has to be noted that SYMPHONY does not include an LP-Solver, but 
can use, through the Osi interface, third-party solvers such as Clp, Cplex, 
Xpress. Furthermore, SYMPHONY also has a structure-specific implementa- 
tions for problems like the traveling salesman problem, vehicle routing prob- 
lem, set partitioning problem, mixed postman problem and others. 

Warm start technology for implementation of Postoptimal Anal- 
ysis. 

Warm restarts are used by modern solvers such as Gurobi, SCIP, CBC, 
SYMPHONY, and others, in order to implement PA. We use the capabilities 
of warm restarts offered by SYMPHONY in order to reduce the solution time 
of the subproblems generated by the computational scheme of LEA [4j . 

SYMPHONY implements warm restarting by using a compact description 
of the search tree at the time the computation is halted. This description 
contains the complete information about the subproblem corresponding to 
each node in the search tree, including the branching decisions that leads to 
the creation of the node, the list of active variables and constraints, and warm 
restart information for the subproblem itself. All information is compactly 
stored using SYMPHONY'S native data structures, which store only the dif- 
ferences between a child and its parent, rather than an explicit description of 
every node. In addition to the tree itself, other relevant information regarding 
the status of the computation is recorded, such as the current bounds and 
best feasible solution. 



3 https: //projects.coin-or.org/SYMPHONY 
* http://www.coin-or.org 
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By using warm restarting, the user can save a warm restart to disk, read 
one from the disk, or restart the computation at any point after modifying 
parameters or the problem data. Note that the use of the PA procedure does 
not always guarantee a positive result. In the case of SYMPHONY, it has 
been observed that, as a rule of thumb, the warm-restart procedure works 
best with a slight change in the conditions of the problem. 



3.3 Benchmarking analysis 

All experimental results were obtained on an Intel Core 2 Duo at 2.66 GHz 
machine with 2 GB main memory, and running Linux, version 2.6.35-24- 
generic. SYMPHONY 5.4.10 was used for the LEA implementation. The re- 
sults are presented in Table [1] where n denotes the number of variables, m 
the number of constraints, k the number of blocks, b the size of separator, 
and underlined is the minimal time of problem solving for appropriate algo- 
rithm. The maximum solving time is denoted by TIMEOUT, and is equal 
to 2 hours. 

Test problem description. All the ILP problems with binary variables 
from a given experiment have artificially generated quasi-block structures. All 
the blocks from a single problem have the same number of variables, and also 
the same number of variables in separators between them. This is required 
in order to evaluate the impact of the PA on the time to solve the problem 
by increasing the number of variables. 

The test problems were generated by specifying the number of variables, 
the number of constraints and the size of the separators between blocks. The 
number and the dimensions of the blocks were calculated by using the number 
of variables and constraints. The objective function and constraint matrix 
coefficients, and the right-hand sides for each of the block were generated by 
using a pseudorandom- number generator. 

Each test problem was solved by using three algorithms, a) the basic MILP 
SYMPHONY solver with the OsiSym interface, b) the LEA in combination 
with SYMPHONY, c) the LEA in combination with SYMPHONY and with 
PA (warm restarts). In all the cases SYMPHONY used preprocessing. 

The computational experiments show that LEA combined with SYM- 
PHONY for solving quasi-block problems with small separators outperforms 
the stand alone SYMPHONY solver (see table [Ij. Additionally, by increasing 
the size of the separators in the problems for the same number of variables 
and block sizes LEA becomes less efficient due to the increased number of 
iteration for solving the block subproblems. LEA's efficiency is improved by 
using warm restarts. The ILP problems corresponding to the same block for 
different values of the separator variables differ only in the right-hand side. 



5 http://coinor.org/download/source/SYMPHONY/ 
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Table 1 Run-time (in minutes) of solving generated problems with quasi-block struc- 
ture. 





Problem parameters 


Solvers 


# 


n 


m 


k 


b 


SYMPHONY 


SYMPHONY -h LEA 


SYMPHONY -1- LEA + PA 


1 


180 


12 


6 


1 


9.503 


0.028 


0.026 


2 


180 


12 


6 


2 


3.019 


0.046 


0.047 


3 


180 


12 


6 


3 


1.671 


0.17 


0.171 


4 


180 


12 


6 


4 


1.164 


0.493 


0.485 


5 


180 


12 


6 


5 


0.084 


5.667 


5.295 


6 


180 


50 


25 


5 


1.572 


0.03 


0.031 


7 


180 


12 


6 


6 


2.7167 


5.321 


5.057 


8 


320 


16 


8 


1 


9.605 


0.025 


0.024 


9 


320 


20 


10 


1 


24.435 


0.029 


0.026 


10 


320 


40 


20 


1 


3.943 


0.016 


0.015 


11 


320 


20 


10 


2 


18.464 


0.092 


0.093 


12 


320 


40 


20 


2 


2.519 


0.048 


0.049 


13 


320 


20 


10 


3 


31 


0.289 


0.289 


14 


320 


40 


20 


3 


3.382 


0.158 


0.157 


15 


320 


20 


10 


4 


47.919 


1.21 


1.204 


16 


320 


20 


10 


5 


14.941 


3.73 


3.682 


17 


320 


12 


6 


6 


44.556 


21.534 


21.636 


18 


320 


20 


10 


6 


23.853 


14.327 


14.043 


19 


500 


50 


25 


1 


90.02 


0.021 


0.022 


20 


500 


50 


25 


2 


TIMEOUT 


0.071 


0.073 


21 


500 


50 


25 


3 


TIMEOUT 


0.681 


0.303 


22 


500 


50 


25 


4 


TIMEOUT 


0.883 


0.879 


23 


500 


25 


12 


5 


TIMEOUT 


7.404 


7.413 


24 


500 


130 


65 


5 


TIMEOUT 


0.076 


0.077 


25 


500 


112 


56 


6 


TIMEOUT 


0.246 


0.244 


26 


800 


120 


60 


1 


TIMEOUT 


0.048 


0.049 


27 


800 


120 


60 


2 


TIMEOUT 


0.119 


0.118 


28 


800 


50 


25 


3 


TIMEOUT 


0.564 


0.563 


29 


800 


50 


25 


4 


TIMEOUT 


2.197 


2.175 


30 


800 


240 


120 


4 


TIMEOUT 


0.036 


0.039 


31 


800 


50 


25 


5 


TIMEOUT 


8.924 


8.638 


32 


800 


50 


25 


6 


TIMEOUT 


31.147 


30.719 


33 


800 


180 


90 


6 


TIMEOUT 


0.399 


0.412 


34 


1000 


50 


25 


1 


TIMEOUT 


0.073 


0.075 


35 


1000 


50 


25 


2 


TIMEOUT 


0.287 


0.286 


36 


1000 


50 


25 


3 


TIMEOUT 


1.07 


1.071 


37 


1000 


50 


25 


4 


TIMEOUT 


3.575 


3.573 


38 


1000 


50 


25 


5 


TIMEOUT 


14.424 


16.327 


39 


1000 


50 


25 


6 


TIMEOUT 


56.359 


59.671 


40 


1000 


250 


125 


6 


TIMEOUT 


0.569 


0.577 


41 


1000 


100 


50 


8 


TIMEOUT 


21.447 


21.414 



These problems can be solved partially by using warm restarts and infor- 
mation obtained from other problems, and this was expected to increase the 
LEA performance. However, the results show a inconsistent behavior. For 
most problems, the warm restarts don't make a difference. For some prob- 
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lems, they improved the solution time (see problems 7, 17, 24), while for 
others they did not (sec problems 5, 15, 40). 

Concluding, there is not a definite answer for the block separator size that 
results to LEA's lower performance (with respect to time to solution) com- 
pared to SYMPHONY. The performance depends on the problem structure; 
as the problem size increases LEA becomes more efficient. 



4 Conclusion 

The main result of this paper is to determine the real computational capa- 
bilities of block elimination algorithms combined with SYMPHONY solver. 
Analysis of the results show that for sufficiently large number of blocks and 
small enough size of separators between the blocks for staircase integer linear 
programming problem the local elimination algorithms in combination with 
a solver for solving subproblems in blocks allow to solve such problems much 
faster than used solver itself for solving the whole problem. 

It seems promising to continue this line of research by studying capabilities 
of postoptimal analysis distributed by other modern solvers. It is also of 
interest to research computational capabilities of local algorithms for solving 
sparse problems of integer linear programming from real applications. 
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